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ABSTRACT 


A theoretical model has been proposed for the problem of a round 

jet in an incompressible cross-flow. The method of matched assymptotic 

expansions has been applied to this problem. Inner and outer expansions 

have been obtained in terms of the expansion parameter e= ^j 0 / r t where 

R. is the jet orifice radius and r is the minimum radius of curvature 
jo t 

of the jet trajectory. For the solution to the flow problem in the 
inner region, the re-entrant wake flow model was used with the re- 
entrant flow representing the fluid entrained by the jet. Higher order 
corrections are obtained in terms of this basic solution. The pertur- 
bation terms in the outer region was found to be a line distribution of 
doublets and sources. The line distribution of sources represents the 
combined effect of the "entrainment" and the "displacement". For 
large values of Uao/U^, the effect of the "entrainment" is small while 
the effect of the "displacement" is important. For small value of 
tWUj, the opposite is true. Pressure coefficient contours for 
Uj/U<» = *25, .5, and .8 were computed for inner and outer region. The 
variation of jet radius, source strength, and mean jet velocity as func- 
tions of distance along the jet trajectory were also computed to aid in 
assessing the role of these quantities in the interaction. 
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I. INTRODUCTION 


One of the methods currently employed to achieve vertical take-off 
performance in VTOL aircraft is to install a lifting jet. An inherent 
problem of such an installation is that it interacts with the crossflow 

O 

due to horizontal flight resulting in an induced flow field which has 
a measurable influence on flight characteristics. 

In recent years, considerable amount of attention has been paid to 
this field. In the experimental area, Jordinson, Ref . [8] made pressure 
measurements in the jet plume which gave jet penetration data and the total 
pressure contours. Later, Bradbury and Wood, Ref. [7], and Vogler, Ref. [9], 
measured the pressure distribution on the orifice plate of the jet. Recently, 
Mosher, Ref. [2], weighed the relative importance of the effect of blockage and 
the effect of entrainment on the induced flow at various ratios of jet exit 
velocity to free cross-flow velocity, i.e., Uj/U ro . The most extensive wind- 
tunnel tests were done by Fricke, Wooler and Ziegler, Ref. [1]. They ob- 
tained data not only for a one jet configuration, but also for two and three 
jet configurations. 

One of the most comprehensive descriptions of the jet in a cross-flow 
was given by Keffer and Baines, Ref. [3]. The picture which emerges from 
their work is of a jet initially comprised of a potential core surrounded by a 
growing turbulent shear region. Within a few diameters of the orifice, the 
jet is completely turbulent, its cross-section becoming kidney-shaped and 
finally developing the appearance of two counter rotating vortices far down- 
stream. Over the initial stages of the jet, the cross-flow is described as 
separating from the jet somewhere near the widest point of its cross-section. 

A "wake" region immediately behind the jet is characterized by negative pres- 
sure coefficients in contrast to the upstream face where pressure coefficients 



become positive. The initial bending of the jet is ascribed to this pressure 
"defect" in the wake. Further along the jet axis beyond the potential core 
stage, the jet is described as fully developed or "established" and some 
tendency toward similarity is noted in the velocity profiles measured by 
Keffer and Baines. The pressure differences in front of and behind the jet 
diminish and the bending of the jet is ascribed mainly to the shear stress 
and momentum exchange due to turbulent mixing or entrainment. 

Among the initial analytical efforts were attempts by Shandorov and 
others. Ref. [12], to derive the jet trajectory. Margason, Ref. [10] has 
proposed an empirically derived equation developed by matching experimental 
data. The results of these investigations have been used by Rubbert, Ref. [11] 
to provide the trajectory at the outset thus reducing the problem to a fixed 
rather than a free boundary one. While Rubbert solves the three-dimensional 
problem completely by computer methods, his potential model contains no 
mechanism to account for the bending of the jet according to the assumed shape. 
Wooler, Ref. [5], actually does derive the jet trajectory, introducing entrain- 
ment parameter and a drag coefficient in the process. These are evaluated by 
fitting his derived trajectory to Margason' s curve. With the trajectory so 
determined, he uses a line distribution of doublets along the trajectory to 
simulate the "local 1 blockage" effect and a line distribution of sources to 
represent the effect of entrainment. This is the so-called sink-doublet 
model. However, this model provides no mechanism to account for the bending 
of the jet and the three-dimensional effect of the trajectory curvature on 
induced flow is neglected. 

Werner and Chang, Ref. [4], solved a complete three-dimensional potential 
flow problem by dividing the flow field into inner and outer regions. Separate 
inner and outer solutions were found in each region and these two solutions 
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were matched by asymptotic expansion techniques in an intermediate region. 

For the inner solutiin, the first order term was found to be the locally two 
dimensional solution to flow about the jet cross-section. For this component, 
the Riabouchinsky wake flow model was adopted. 

For the outer solution, the first perturbation was found to consist not 
only of the familiar doublet distribution but a source distribution as well. 

This latter distribution resulted from the displacement effect of the jet 
as turned in the downstream direction. As such, it had two components, one 
proportional to the trajectory curvature and a second contribution due to 
the increase in jet cross-sectional area. This result raised a rather funda- 
mental question about the use of net sink distributions to represent entrain- 
ment effects: If the fluid entering the "sink" is to be considered a part 

of the jet, the jet displacement increases requiring an additional source 
term in the solution and the net sink strength should be zero. 

The flow model of Werner and Chang, Ref. [4], did not consider explicitly 
the effect of entrainment by the jet. It's the purpose of this investigation to 
explore the analytical structure of a model in which the effects of entrainment, 
blockage and curvature exist concurrently. The particular model has been 
chosen largely an mathematical grounds and admittedly an element of artifi- 
ciality has been introduced thereby. On the other, hand, it provides a 
framework in which to explore the mutual dependence of entrainment, trajectory 
curvature and blockage, while the resulting insight should be useful in the 
treatment of more rigorous models. 
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II. THE FLOW MODEL 


« 


The model to be used in this investigation together with pertinent 
coordinate systems are illustrated in Fig. 1 and Fig. 2. The external 
cross-flow is assumed to be steady and inviscid. A vortex sheet which 
encloses the whole jet is the inner boundary of this problem. The plane 
normal to the jet trajectory intersects the vortex sheet to give a shape 
shown in Fig. 4. For external cross-flow, this vortex sheet can be re- 
garded as a streamline in the normal plane. The streamline (or vortex line) 
CABHF is assumed to separated from the jet at B, and a wake of constant pres- 
sure Pw(s 1 ) is formed immediately behind the separation point. Then this 
streamline rolls-up, forms a backwards stream HH / FF / at the center of the 
wake and flows into the jet. The backwards stream HH , FF / is called re- 
entrant flow. 

The geometrical shape ocf the wake and flow-rate of the reentrant flow 
is governed by the cavity number Q and Uoo, where 


PW - Pa, 



(II-l) 



U = U_ sin a 
n 00 


(H-2) 


If Q is assumed constant along the jet trajectory, the shape in each 
normal plane will be geometrically similar. In this investigation, we do 
assume that Q is constant. 

In each normal plane, the reentrant flow carries some fluids of external 

I 

cross-flow and flows into the jet. Then, these fluids are carried away by 
the jet into the next normal plane and forms so called entrainment. Two 
assumptions have been implied here; all entrainment occurs in the wake and 
there will be no entrainment if the external cross-flow velocity is zero. 
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But from the study of Keffer and Bains (Ref. [3]), the entrainment of a 
jet in a cross-flow is much greater than that of a free jet. And, the 
backwards flow in the center of wake of jet is confirmed by Mosher 
(Ref. [2]). These assumptions might not be bad after all. 


III. MATHEMATICAL FORMULATION 


a) Formulation of the problem 

Since the flow exterior to the vortex sheet (Fig. 1) is assumed to be 
inviscid and steady, a velocity potential $ exists and is governed by La- 
place's Equation. A boundary valued problem is thus formed as following: 

V 8 § = 0 (III-l) 

1. Inner Boundary Condition: On the vortex sheet the velocity 

vector is tangential to this sheet. 

2. Outer Boundary Condition: At infinitely far away from the 

jet, the velocity vector is i . 

3. On the plane Z 10 = 0, the velocity vector is tangential to this 
plane. 

b) Image System 

The third boundary condition in preceding section is solved by re- 
flecting the flow field into the lower half space (dotted line in Fig. 2). 

Now, the new flow field is symmetrical with respect to the plane Z lo =0; 
the velocity on this plane is, therefore, tangential to this plane. 

c) Inner and Outer Solution 

The method of asymptotic expansion (Refs. [13, 14 and 15]) is used 
to solve the problem formulated in preceding sections. By this method, 
separated "inner" and "outer" solutions are introduced, with inner solution 
satisfying the inner boundary condition and outer solution satisfying the 
outer solution. In a intermediate region, these two solutions are matched 
asymptotically. Also, Laplace's Equation is satisfied by both inner and outer 
solutions. The use of this method is prompted by the fact that, due to three- 
dimensional nature of the problem and geomtrical complexity of the inner boun- 
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dary, it is extremely difficult to obtain an analytic solution which can be 
applied throughout the flow region. And also, for iWUj « 1, the radius of 
curvature r fc of the jet will be large compared to the jet radius Rj, i.e., 
Rj Q /r t << 1. Under this condition, the flow near the jet, within a radius 
of the order of Rj, is expected to be nearly two-dimensional, but not com- 
pletly so. 

Introducing 

R. 

e = (III-2) 

t 

as a parameter, we write inner and outer solutions as asymptotic series: 
i 00 i 

cp (X,Y,s; e )~£ g m (e)cp m (X,Y,s) (III-3a) 

with X,Y fixed as e -» o 

CO 

cp°(x,y,s; e)~ £ f m ( £ ) <P° ( x > Y> s ) (III-3b) 


with x,y fixed as e-*o 

where 

$ 

cp = — non-dimensional potential (III-4a) 

*-'00^ £ 


X = 





non-dimensional inner variables 


(III-4b) 


x 




1 

s= non-dimensional outer variables 

r 

t 


(III-4c) 


f 's and g 's are functions of the parameter e only with the property: 
m m 


£ m (e) 


-• O 


as 


G - 


(III-5) 


7 


i 


Eq. (III-4) establishes R as the basic normal length scale and as 
the basic tangential (to the trajectory) scale for phenomenona in the inner 
region. The basic length scale in the outer region is r for every direc- 
tion. At this point, attention should be called to the following fact: 

3§ 


••••]- •[*£] 


~ u “ [ f o (e) 3 s° + ] °[ f o (e) ] 




(III-6a) 


(III-6b) 


If f Q (e) is taken to be unit, which will not lose generality in comparing 
order of magnitude, then 3$/3x i is of order (l/e)> while 3$/3s i is of order 
(1). In other words, following important assumption is implied in this way 
of formulating inner and outer expansion; in inner region, variation in $ 
for a displacement of o(R^ q ) in X-Y plane are of the same order as variation 
for a displacement of o(r t ) in s direction. This assumption is, of course, 
true as long as ILo/U^ remains small, 
d) Basic equation in dimensionless variables 

As stated before, both inner and outer solution should satisfy Laplace 
equation. In inner region, Laplace's equation in inner variable takes the 
form (Ref. [4]). 

-^cp 1 + 5 s cp 1 _ e h(s) 3c^ sfScp, 1 JL( I > + 

3 X 2 3 Y 2 1 -ehx 3x e L5s 3s V 2(1 - e h X ) 2 J 

( i y 3 2 cp i i = o 

\1 - e hx/ 3s 2 J U 


(HI-7) 


where 


h(s) = r t /r 1 (s 1 ) 


Expanding the coefficients into power series, Equ. (III-7) can be written 
as 
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Ix^ + T& ■ eh(s) %x~ + £S [t^ ■ hS(s)x lx . 


+ 0(e 3 )~ 0 (III-8) 


Since the specific form of S m (e)' s *- n E< 1* (III-3a) is still unknown, we 
cannot substitute Eq. (III-3a) into (III-8) to obtain the governing equa- 
tion for each cp 1 's. 

In outer region, the parameter e does not appear explicitly in the 
Laplace Equation in outer variables. Thus, Equ. (III-3b) can be substi- 
tuted into V 2 cp = 0 to obtain the governing equation for each cp° as: 


V s cp — 0 , for m ^ o 

mb 


(IH-9) 


where V 2 is the three-dimensional Laplacian. 
e) Boundary conditions 

If the vortex sheet (Fig. 1) is given by the equation, 

F(x x > y 1 , Sj_ ) = 0 

the over-all inner boundary condition will take the form: 

Vcp 1 • VF = 0 

But the inner boundary condition for each cp* is yet unknown. 
The over-all outer boundary condition is 


(III-10) 


(III-ll) 


lim „ o — 
Vcp = 1 
x,y- » ^ 


(III-12) 


By substituting Eq. (III-3b) into this equation, the outer boundary con- 


dition for each cp° are obtained 


o t 


lim „ _ . 

Vcp = i 
X,y — a> ^ o 


(III-13a) 


lim 


Vcp = 0 , for m s 1 
x,y-*a> ^ m 


(III-13b) 
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The remaining conditions are supplied by requiring the inner and 
outer solution to match asymptotically over an intermediate region. 
Formally, this can be stated as 


lim 
e -• o 


[" -[“ £ n <e) »n] 


0 


y(e) 


(III-14) 


The notation [ 1 denotes behavior near the outer bound of the inner 

region and [ ] ^ the behavior near the inner bound of the outer region. 

Eq. (III-14) states that over the match region, the M*"* 1 partial sum of 
cp* agrees to the N*"' 1 partial sum of cp° to the order of [y(e)]. 
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IV. SOLUTIONS FOR cp° AND m 1 

As e-»o, the disturbance vanishes altogether, the outer flow re- 
gion is just the undisturbed flow, which is independent of e. These 
conditions will be satisfied if we take 


f (e) = 1 
o 

f (e) -* o, as g -* o for m 2 1 


cp Q (x,y, s) = 


cosq(t) dT + xsina(s) 


o 

= co° (o,o,s) +x sina(s) 
T o 


(IV-1) 

(IV-2) 


(IV-3) 


where cp°(x,y, s) is the undisturbed potential in curvilinear coordinate, 
o 

The matching condition, Eq. ( III -14) for m = o, n = o becoming 


lim 
e -• o 



cp° (o,o,s) + eXsina(s) 
Y(e) 


0 


(IV -4) 


It is obvious that we should take y(e) = 1> then this gives 
g o (e)=l 

s 

fcp 1 "! =cp°( 0 , 0 ,s)= r COSa(T)dT 
L J o o 

Now, the inner solution stands: 

i i . , x i , 

cp ~ cp Q + gi(e)cpi + ... 


(IV-5) 


(IV-6) 


Substituting this equation into the basic differential equation (III-8) 
and basic boundary condition (III-ll) gives 


*o " ° . 

/dcp v 

7 xy< ' V " \—> * 0 


(IV-7a) 
(IV- 7b) 
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where and are two-dimensional Laplacian and two-dimensional 
gradient operator respectively. Together with the outer boundary condition 
given by (IV-5a), the governing equations for cp^ are completely deter- 
mined. The solution of this set of equations is simply: 

CPp = cp° (0,0,8) (IV-8) 
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V. Solutions for cpi, Cpx and cps 


a) General solution for cpx 


cp x has to satisfy following set of equations: 


V 3 cp?= 0 


lim cpi = 0 
x,y-*» 


(V-l) 

(V-2) 


With the inner boundary condition along the jet trajectory unknown. By 
potential theory (Ref. [19]), the most general solution satisfying these 
conditions can be represented as a infinite series of a line-distribution 
of singularities along the jet trajectory: 


„ dT 7 " (1> 

CP! = I ^ + | 


0 

-CO 


-|a; y (T) dT ” \ (T ) COS y(T ) dT 

— + 1 -s 

-co O 


I 


I- 1 a < T ) p 2 [cos y(t)] dT 


• • • • 


(V-3) 


The first term is a line distribution of sources, the second term a line 
distribution of doublets, the third term a line distribution of quadrupoles 
and so on, with the symbols having following meaning: 


^ 1} (T) 




= T|(T), C(T)] Strength of source 

= Ui^C^CO, ri (T ) , C(t)] Strength of doublet 

= ni^C^CT), n(T), C(T)] Strength of quadrupole 


r „ - 1 < x o-?> + j + k cz -c> 


r = r 
o * o 


cos y (s) = i • — 


(V-4) 
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X o = x * /r t’ y o = y i /r t’ Zl “ Z 1 /r t 

x , y , Z : location of the point where the potential is to be evaluated 

o o o 

?» T)» C : location of the singularity on the jet trajectory 

^n : nt ^ or< * er l*®rgrangian polynomial 

The strength and so on will be determined by matching 

with inner solutions. 

In order to do the matching, the behavior of cp°, near inner boundary 
has to be determined first. In other words, we need the asymptotic expan- 
sions (with X, Y fixed as e -* o) of each singularity (source, doublet, 
quadrupole and so on) . 

Fortunately, only expansions of sources and doublets are needed for 

our work. These expansions are obtained from Werner and Chang (Ref. [4]). 

In all following formula, the index i means the expansion is obtained from 
til 

the i tn order outer solution. For example, in this section we are dealing 
cp°, then we should take i = l. See Fig. 3 for symbols used. 

1) Asymptotic expansion of a line distribution of sources 


where 


lim ^ -H^COdT 

X, Y fixed j 

G — ♦ O -oo O 

~(loge)2n^ ) (s) +[2|j< i) (s)log (X 2 + Y 2 )^ +T.(s)] 
+ (e log e) ^^(s) h(s) X 

+ cx[nj l) (s)h(s) + ^ l) h (s) log (X s +Y 2 )^+P i (s)] 
+ 0 (e 2 log e) 


T^s) 


=[J + J] 

- co A 


•|Jl^(T)dT 

C 


- 2^ i) (s) log 2 6 


(V-5) 


(V-5a) 
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(V-5b) 


- 6 


— f» p 

P.(s) = + 

x L w v _ 


-H^CO cos p (t) dT 


-U^ l) (s)h (s) log 2 6 


6: a small number, which has the property: 6 -» o, as e -* o 

C: the chord of the trajectory between points P andP' (Fig. 3) 

P: the angle between chord c and x-direction (Fig. 3) 

This expansion is correct only if T^(s) and P^(s) are independent of 6. 
But this can be easily proved. 

2) Asymptotic expansion of a line distribution of doublets 


lim 

X, Y fixed 
e -♦ o 


CO 

r 


co 


[i 1 1 ^ (T ) COS y(T) dT 


(V-6) 


, 2|ai 1 \s) cosQ m 

~ ~ T (log ; (s) h (s) 

e (X 2 +Y 2 ) 2 

+ £i 1 (i) (s)h (s)^|ie _ u ( i > (s)h (s) log (x 2 + y 3 )^ 

+ G. (s) 

+ o(e) 

where 

cos 8 = — x 

(X s + Y 2 ) 2 

- 6 <*> (i) 

G.(s) =[ J + Jj -p-j (t - s) 2 cos(a- a ') 

- CO 6 

-Kt - a) sin (a- a' jdT +,u ^ 1 \s)h (s) [log 2 6-|] 
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(b) Governing equations for cp i and cp\ 


The inner and outer expansion stands 


O O rr -^o l)(T)dT p U ( i 1 ) (t)cos y dT 

cp ~co o (o,o,s) + X sina+f 1 (e)!J + J -ps 

o o 

— 03 — 00 

M.2 1 ' > P 2 Ccosy> 

+ j _ + J + ... 

O - 1 


with x,y finite as e-*o 


(V-7) 


cp 1 ~ cp°(o,o,s) +g x (e)co\ + . . . . 

with X,Y finite as e -* o 


(V- 8) 


In matching region, we have 


£ f n/ e ^ ^m = ‘ 0 °’°’ s > + eXsin0! + f i (e) ( lo ge) 2 ^/ 


1 f 2 ^ \ 1 /' 2kl ® 1 \ 

+ 7 Re \ Z / + e" 5 Re l” ) + — 


( 1 ) 


(V-9) 


[sS m (e)c Pm] = c Po (o, ° >s)+ 8 l(e) [ c Pi] 


(V-10) 


For later convenience, we use here the complex variables Z and z: 


Z = X + i Y 
z = x + i y 


Re(Z ) = X 
Re(z) = x 


lm(Z) = Y (V-lla) 

Bn(z) = y (V-llb) 


By taking y(c ) = f x (e ) and being simplified by the assumption e/gi(e)-*o as 
e -* o (we will varify this assumption later), the matching condition 
Eq . ( III- 14) becomes: 
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,( 1 ) 


e* <£> [< l °S e > 2 ^ 1)+ r 

- Si (e^cpi 1 ] } = 0 


2H 


<1) 


which gives: 


g i (e) = f i (e) log e 

[cpi] = Zii^^Cs) 


o 

( 1 ) . 


= 0 , f or P > 1 




(V-12) 


(V-13) 
(V - 14 ) 

(V-15) 


Although (V-14) gives the outer boundary condition for cpi, the specific form 
of g x (e) is still unknown, and therefore we are unable to obtain the governing 
differential equation for cp j . A lot of information is revealed by (V-15): 
it says the strength of all higher order poles are zero and the first, order 
outer solution is a line distribution of sources only: 



(V-16) 


This is just the device used by many investigators (Ref .[5,6]) to simulate the 
effect of entrainment. They generally used an empirical formula for the 
strength of the sources. In this investigation, the strength of sources, 
which is not yet known* 'will be derived as an integral part of the whole 
theory. 

We proceed to next order. In the matching region: 


rj f m (e) cp° J = cp° (o,o,s) + e X sina + f x (e;)[d°g e) + 

+ 2n ( o 1) (s) Re (log Z) + Tj (s)] 

Again, using matching condition (III-14) with y(e)=g 2 (e), we have 


(V-17) 
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(V-18) 


Ss (e) = f i(e) = e 

[cpa"l = x sin CH- 2 |_l^^ (s) Re (logZ) +T X (s) 
g i (er) = g 3 (e ) log e = e log e 

Previously, it was assumed that e/g i (e ) -* o as e -* o . Looking at (V-20), 
we see that this indeed is the case. 

Now, the inner solution stands: 


(V-19) 

(V-20) 


cp 


1 i , , i , i , 

~ cp + e log ecp! + e cp 2 + . . . . 


(V-21) 


Substituting this into the basic differential equation (III-ll) gives 


v^cpi = o 




i \ 

7 xy^ ■ V ■ b 7) a -° 

Sep 1 

7 xv t( ’ s ■ h(s) -ax' 

i /SepJ . 

V'f’s • 7 XY F ■ vsr) ■ 0 


(V-22) 

(V-23) 

(V-24) 

(V-25) 


Here, cr=a.+a is inner boundary, see Fig. 4. 

J w 

c) Solution for cpi 

From Eqs. (V-22), (V-23), and (V-14), cp\ is governed by: 




= 0 


= 0 


(V-26a) 


(V-26b) 


lim 


X,Y — < I , » " 2 ^o 1)(8) 


(V- 26c) 
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The solution of this set of equations is simply: 

cp i (X, Y,s) = 2^ 1) (s) (V-27) 

Therefore, cp^ represents a flow field which is parallel to the jet tra- 
jectory. At each station s, the velocity of this flow field is constant 
throughout the whole inner region, 
d) Solution for cpg : The Re-entrant Flow Model 

Since cp^ = cp°(o,o,s) is independent of X, Eq. (V-24) reduces to ^j^cpg = 0. 
Together with (V- 25) and (V-19), cpg is governed by following set of equa- 
tions : 

VyyCpa = 0 (V-28a) 


x y^cp s = X sin a (s) + 2|j^ 1 ^ (s) Re (log Z) 



(V-28b) 

(V-28c) 


Here, |j^^(s) is still unknown and will be supplied by cp| . 

A number of two dimensional wake models can be chosen here to give the 
solution of this set of equations, each will give a different inner boun- 
dary shape. The one we choose in this investigation is the "Re-entrant 
Flow Model" (Refs. [16,17,18]). The flow field of "Re-entrant Flow Model" 
is illustrated in Fig. 4. The stream line CABHF is separated from the jet 

at point B. B'AB is denoted by a. and B’H'F, BHF by o . a - a. + a con- 
r j wjw 

stitutes the inner boundary. The shape of the rear surface of the jet 
(dotted line BFF'B') has no direct effect on the external flow since the 
external flow is bounded by a rather than by the rear surface of the jet. 
The front surface of the jet (a.) is assumed to be a circular arc. The 
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geometry shape a w and the wake pressure Pw are characterized by cavity number 

Q. The actual shape of a = CT. + a for Q = -1.5, -2.0 and -3.0 are demon- 

J w 

strated in Figs. 6a, 6b and 6c, respectively. 

Mathematically, the flow region of the "Re-entrant Flow Model" consists 
of "two sheets of Reimann surface". What we have just described is the first 
sheet. The re-entrant flow HH 'FF' carried some fluids of the first sheet flows 
into the so-called second sheet. The flow field of the second sheet is shown 
by the dotted lines FE and F'E'. It should be noted that the second Reimann 
sheet is required only for mathematical consistancy of a two-dimensional wake 
flow problem. In context of our three dimensional theory, all fluids which 
flow across line FF' (Fig. 4) will be carried away by the jet to next plane 
further along the jet, where this fluid will cause an increase in cross- 
sectional area of the jet. Thus, in present formulation, entrainment is re- 
presented by the re-entrant flow. 

e) W a (Z,s): The Solution of Re-entrant Flow Model 

Solution to the "Re-entrant Flow Model" has been well documented 
(Ref. [16, 17,18]). This solution has been obtained by the method of Levi- 
Civita for cavity flows about curved surfaces. Only a summary is presented 
here. In two dimensional problem, it is convenient to work with complex 
variables. We are looking for the complex potential: 

W 2 (Z,s) = cpa (Z,s) + iy a (Z,s) (V-29) 

Here Y 2 is the stream function. By using conformed mapping technique, both 
W 2 - plane and Z - plane are mapped conformed onto an auxilliary plane; the 
Q - plane. The shapes of the boundary a = cr^ + a„ in Z-plane, w 2 -plane and 
£-plane are illustrated in Fig. 5. The boundary a in the Z-plane is trans- 
formed into a geometrically simpler shape in the £-plane, a semi-circle. The 
solution can be represented in a parameter form as: 
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1 dws / 

TliTa ~ ( ^ } = 7Ina’ (U2 ~ iV2 


1 - -1V„)-(1-Q)* (C-id)(dC- l ^ . -i) e F(C) 

u Q; (C+id) (dc+i) (C+i) 


R. s 
J 


dw 

1 8 (r)-b(l-o^ (S 4 -D (C g +d 8 ) (d 2 c 3 + l) 

iina dr u C(C 2 +c 2 ) s (c 2 r 2 +l) 2 


(V-30) 


(V-31) 


w . 


R. sin a 
J 


(C) 2 c 4 (1 -Q>* lo s[ (£2+c 2 )(C 2 + l/c) *] 


_ b n _ 0 )l (c 2 -d g )(l-c 2 d 2 ) C 3 + B 

bU Q) 2c 4 (£ 2 +c 2 )(£ 2 + 1/c 2 ) o 


where 


B o=i 




k (c s -d 2 )(l - c 2 d 2 ) 
2c 4 (1-c 2 ) 2 


(V-32) 


(V-32a) 


F(£) is Levi-Civita's Function: F(C ) = iA a Q - iA 3 Q 3 /3 + iA 5 £ 5 /5 + . . . 

Here, the constants A 1 ,A S ... will determine the shape of the line AB (Fig. 4). 
For this investigation, only two terms of F(£) are taken to approximate the 
Levi-Civita's Function, 


F(C) = iA lC -iA 3C 3/3 


(V-33) 


To obtain Z as a function of £, w s is eliminated between (V-31) and (V-30) 

_L Ob (r) - dWs/dC =b (C+id) 3 (d C+i) 3 (C+i) (C 8 -1) e [-iAxC+iA3C 3 /33 
R, dw 2 /dZ C (C S + c 8 ) 2 ( c 3 C 2 + 1/ 


(V-34) 


All together, we have bjCjdjAjjAg five unknown constants. If these five 
constants are known eq. (V-34) can be integrated numerically by choosing a 
path in the Q -plane and integrating step by step to obtain corresponding 
values of Z. It is desirable to extend the integration to values of z 
well into the matching zone which for small values of Uoo/U^ may as touch as 
five to ten diameters from the origin. Unfortunately, for this range of 
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values of 2 , Q must approach the neighborhood of the singularity c in 
the £-plane (Fig. 5), where small incremented in Q yield unmanageably large 
ones in Z . To avoid this inconvenience the reverse procedure of first 
choosing a point in the Z -plane and determining the corresponding values 
of £ by a modified Bairstow scheme (Ref. [20]) was employed. 

The constants A x ,A 3 ,b,c,d are determined by following five conditions 
(Ref. 16): 

1) a predetermined cavity number Q 

2) a predetermined curvature at A (Fig. 5) 

3) the curvature at B must be finite (Fig. 5) 

4) the curvature at B must be equal to curvature at A 

5) Z(£) must be a single valued function at Q = ic 

After applying above five conditions, the following five equations result 


(c+d)(cd+l)(c+l) (A i c+ A 0 c/3) 
(c-d)(cd-i)(c-i) -u-y; 


2 + 2d 


2d 


1-c 3 ' cd+1 " c(c+d) +A i +A s cS 0 


(V-35a) 

(V-35b) 


x£p+ 1 + =° 


(Af 3A s )(l+c 2 ) 4 (Aj + A, )(l-c 8 ) 4 


4 ( 1+d 2 ) 


2 \2 


8 ( 1+d ) 4 


|“ A i ~ a 3 /3) 


( A 1 +A 3 )(1-c 8 ) 4 (_a 1 -A 3 /3) 
8 (1+h) 4 6 


(V-35c) 


(V-35d) 

(V-35e) 


Of course, this is another set of nonlinear equations. By suitably constructed 
numerical methods, this set of equations is solved to obtain the five constants 
for a predetermined cavity number Q. A x and A 3 determine the shape of CT^. 
b is a scale factor', ic and id are the points in Q -plane corresponding to 
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infinity and rear stagnation point in Z -plane respectively (see Fig. 5) 


The procedure outlined here gives only an approximated circular arc along 
o^. The constants A x and A 3 are chosen so that the curvature at front stag- 
nations point A and points of separation B and B' will have the same value. 
From actual numerical computation, the maximum variation of curvature be- 
tween A and B is less than 2% for Q=-2.0. i n viewing the extreme uncertainty 
of the shape of the jet, this is more than enough. In Fig. 6a, 6b and 6c, the 
actual shape of a for Q=-1.5, -2.0 and-3.0. 


f ) Power Series Expansion of w 2 and the strength of sources jj ^ 
w 2 can be expanded into a power series having the form 

a (s) 

w 2 (Z, s) = sina(s) Z + b(s) + a (s') log Z + £ — (V-36) 

OO ^ £1 

It should be noted that Re[w 2 (Z, s) - b Q + T x (s)], not Re[w 2 (Z,s)], is the 
solution of the set of eqs. (V-28). From this, we obtain the behavior of 
cp 2 at outer boundary of the inner region as: 

[cpal = J" Re[w 2 (Z,s) - b Q (s) + T X (s) 

L J o L J o 

= X sinCt(s) + T X (s) + a Q (s) Re (log Z) (V-37) 


Comparing (V-37) and (V-28b) gives 

(1) 


a (s) 

^o (s) = “V" 


(V-38) 


Thus, the strength of the line-distribution of sources for outer expansion 
is determined provided that a Q (s) can be evaluated. 

The evaluation of the constants b Q , a Q , ... in (V-36) presents 
some difficulty because w a (Z,s) is not an explicit function of Z. At least 
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three constants b , a , a, have to be evaluated fort/ 1 is) and for the later 
o o " o 

matching with outer solution. The rest of this section is devoted to obtaining 

these three constants. 

i) a (s) 
o 

Differentiating (V-36) with respect to Z gives 


dw . 


a a • 


dZ sina + z - z3 + .... 


(V-39) 


By Residu theorem 


. S o _ f 1 dw „ _ f 1 

17 L R . s in a jR.sinadZ j R.s 

J r' J •p J 


dw 

in a d Q d C 


(V-40) 


where r is any circle which encloses the surface a in Z-plane and Ti is the 
image of F in £-plane (Fig. 5). After expanding dw 2 /d£ (V-31) into power 
series, (V-40) gives 


= . b hf 

R. sin a 2 c 4 ' 

J 

ii) a^s) 

Multiplying (V-39) by Z and using Residu theorem again, we have 

2 . ai l dWg z dz r l dWs z l dz 

m R?sinQ J sin a dZ R. R. J sin a dZ R . R . df 


(V-41) 


J J 


J J 


(V- 42) 


Next, expanding dw 2 /dZand dZ/d Q into power series and integrate dZ/d£ term 
by term: 


sin 


i dw 2 

iTa dz” = b 0 + b i ( £ " ic)+ '** » for |c- ic |<U-c) 


(V-43) 


<c) - ( g T ic)" - £ c i + C2 ( c - ic > +c 3 (c- ic) 3 + ] 

for 0 < j Q- ic | < c 


(V-44) 
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(V-45) 


R <C> ~ q . ic + C Q + c 3 <C‘ ic > + ••• 

3 b 

for 0 < |£- ic | < c 

where c is a constant of integration and c 2 = o since Z/R^ is a single 
valued function of Q. Substituting (V-43), (V-44) and (V-45) into (V-42) 
and simplifying it, we have 


a i 

n o . = b , c Ci 

Rrsma 1 o 1 
J 


- b. 


(V -46 ) 


where 


(c+d) a (cd + l) g (l+c) s (1 + c 2 ) (Aj c + A 3 c 3 /3) 
1 4 c 3 (1 -c 4 ) 2 6 


(V-46a) 


c 

c = Um f 
O ^ ic J 


1 r, (c+id) a (dC + j) s (C+i) 2 (c s - 1) (-A,C+iA,C3/3) 
(Q- ic ) 2 L C(C + ic > 2 (c 3 C 2 + l ) 2 


_ ic x 

C iJ d £ -<— > - 1 


(V-46b) 


b 



+ A 3 c a 



e (A x c-A 3 c 3 /3) 


(V-46c) 


b 2 = (1- Q)*{fi [2A3C+CAJ +A 3 c 3 ) 2 ] +2 if a (A x +A 3 c a ) + f 3 } 

~ ("A x c-A 3 c 3 /3) 


f - (.e ; .d )(cd-l)(c- _ l) 
1 1 ~ (c + d)(cd + l)(c + l) 


(V-46d) 

(V-46e) 


_ ~ 2 i r d(l+d g ) d(l + d 2 ) (1 + d) 5 1 

f 2 (1 - d) a L (c + d) 2 (1 + cd) 2 “ (1 + c) 2 J 


(V-46f ) 


4 r d(l+d g ) . d 2 (1 + d 2 ) (1 +d) 2 ~| 

3 ~ (1 - d) 2 L (c + d) 3 (1 + cd) 3 (l + c) 3 J 


(V-46g) 


Ihe evaluation of c Q (V-46b) deserves some explanation. By suitable numerical 
method, (V-46b) is evaluated to give a sequence of numbers as Q approaches ic. 
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This process continues until the sequence converges to a stationary value, 
which is just c q . Actually, if Q is too close to ic, the sequence will 
diverge eventually because of round-off error of the machine. Neverthe- 
less, a stationary value is clearly observed before its final divergence. 


iii) b o 
Examining 


b 

o 


R . sin a 
J 


(V-36), we can obtain b Q by the limiting process 


lim 


nT w 2 

» UTTlina < z) - 


R . R . s in 
J J 


a los Z ~ 


Q- ic^R. sina ^ ~R. ^ " R. sina ^ log R. + 1 ° s R j)] 

5 J 1 ■* 1 J 


(V-47) 


For |£ - ic| < (1 - c) and ]£ - ic | < c, w 2 (V-32) can be expanded into a power 
series having the form 


w a 

R. sina ^ 
J 


<*i 

(C - ic ) 


+ d 2 log(£ - ic) + d 3 +d 4 (£ - ic) + .... 


(V-48) 


It is easy to show that 

d = _ c = - b(l -0)^ ( cg - d 8 )U - c 2 d 2 ) 
di cj bU q) 4c 3 ^ _ c 4) 


a ^ 

^ R. sma 
J 


A - X- 

u 3 R-. sin a 




b (c 2 - d g )(l - c 2 d 2 )(3c 4 +8c 2 +l) 


8c 4 (l - c 4 ) 2 


(V-48a) 

(V-48b) 

(V-48c) 


Inserting (V-48) and (V-45) into (V-47) and after some algebraic work, we 
arrive at: 


b a R. r / 1 s \ s i R. 

2 = —£ 12 loe T ~ cS ) 3 L _J2 1 

Rj sina sina R^ S [_2(c 2 - d 2 ) (1 - d 2 c 2 ) + R^J 

(V-49) 

+ b (l -0)^ (c s -d 8 )(l-c s d 2 )(3c 4 + 8c 2 -H) . in io_ ^jo 
DU ^ 8c 4 (1 -c 4 ) 8 lTT sinaR j 
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g) Entrainment E g (s) and drag coefficient (s) 

The entrainemt E 2 (s), flow-rate per unit length of the jet, is assumed 
to be the two-dimensional flow-rate of the re-entrant flow, which is just 
the sink strength -a^ multiplied by 2tt : 

E 2 (s) — a d s i 

Rj(s)sina(s) Rj(s)sina (s) TT ^'^3(1“Q) (V-50) 

Since b,c and d are functions of Q only, E/R. sin a is also a function of Q 
only. 

The drag D 2 (s) on the jet per unit length can be easily obtained by 
considering momentum balance of a two-dimensional control vol um e which 
encloses the normal cross-section of the jet (Ref. [16]); 


D 2 (s) = p E s (s) (sina(s) + U^) 


(V-51) 


and the drag coefficient: 


Eg (s) 


C D 2 i p 14? 2 R. 


Rj (s) sin a (s 


) ( 1 + sina ) 


(V-S2) 


where U w is speed of the free streamline BHF(Fig. 4); U w = (1 - Q) 2 sin a . 
C D gis also a function of Q only. The graph of both C^ 2 and Eg/R^sina are 
presented in Fig. 7. 
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VI. SOLUTIONS FOR cp ' % , cp 3 , cp 3 , cpi and cp 5 


. o ,0 

a ) CPs and cp 3 


Now, the inner and outer solutions stand 

“ -y(T)dT 

cp° = cp°(o,o,s) + x sina(s> + g | + f s (e) CP2 +f 3 (e)cp 3 + . .. 


with x,y finite as e -> o 


(VI-1) 


os a 


cp 


1 =cp°(o,o, s) + (e log e) a (s) + e R Z sina(s) +a log Z+T 1 (s)+v; 
~ o o e l o 


1. i i 

+ g 3 (e)cp 3 +g 4 (e)cp 4 + g 5 (e)cp B + ... 


with X,Y finite as e -* o 


(VI-2) 


As stated before, the most general solution for cp 3 can be written as: 

^ -|ag(T)dT Tu i (T)cos y(x)dT Z p \ (t )P 3 (cos y)dT 

cp° -J + J Js + I p 


+ . . . 


(VI-3) 


After inserting (VI-3) into (VI-1) and obtaining behaviors of inner and outer 


solutions at matching region, we can have the matching condition (III-14) as 

«.te)T l ” 


iim 
e -* o 


f 2 (e ) (log e) 2 Pp Re (-^ 7- L ) + 7 * Re (~Z^ 2 ) + •** 


y(e) 


cp 3 


fl = 0 (VI-4) 


By taking y(e) = f 3 (e) log e, we have 


g 3 (e) log e 

[q>i] o =2u»( s ) 


(VI-5) 

(VI-6) 


U 8 = 0, for P s 1 
P 


(VI-7) 
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Just as cp° , cp ° is also a line distribution source only. Actually, cp° is a 

correction to cp° . In Chapter V, we obtained the local velocity at infinity 

in each plane normal to the jet trajectory to be sin a and this gives 

on as a line distribution of sources with strength a (s)/2. We will see 
T o 

later that the local velocity at infinity sin a is only correct to the order 
of e in inner region. If we go to a higher order, the local velocity at in- 
finity has to be modified and this will, in turn, cause a correction cp^ to 
o . 

cp! in outer solution. 

We proceed to cp° with the same technique; this time we have: 

f 3 (e ) = e 2 

a j (s) 

H|(s) 2 


(VI- 8 ) 

(VI-9) 


|i|(s) =0 , for Ps 2 


(V Ir 10) 


(VI-9) and (VI-10) gives 


00 



- 00 


(t )dT 


r 

o 


+ 


00 sLi 

P 


oo 


(t)cos y (t ) dT 



o 


(VI-11) 


Thus, C 0 ° is a line-distribution of doublets and sources with the strength of 

sources |j3 unknown up to know. The line distribution of doublets here is 
o 

just the device used by many investigators to simulate the effect of blockage. 
Later, we will see that the source term represents the three-dimensional effect 
of curvature and growth of the cross-section of the jet. 


b) The Governing Equations for cp^ , cpt anc * cp 5 

In the matching, the inner and outer solution (VI-1) and (VI-2), with the 
aid of (V-5) and (V- 6 ), become: 
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p f m (e)cp° J ~cp°(o,o,s) + e X sina+ ep Q loge +a Q Re (log Z) + T X 

a , a a , _ 

+ (e logo) -fhX + eX^y h+-^ hRe (log Z ) +P x J + 0( e 8 )J 


+ f 2 (e)plog e) 2|a® + 2|j 2 Re (logZ) + T 2 +0(e loge)~j 

+ e 2p a x Re(i)+aoge)(2 U 3 -^h) + ( 2 ^ - h ) Re (log Z ) 

+ T 3 +G 3 +-§- (| + |) +0 <e log.e)] 
p g m (e) ~ cp° (o,o, s) + (e log e) a Q +e Rep- sina + a Q log^ + T! 


(VI-12) 


+ S ^“Tn 2 '! +[f 2 (e) log e] 2|a s +g 4 ( e )["cp4 +g s (e)[ c P5l (VI-13) 

Inserting (VI-12) and (VI-13) into the matching condition (IV-14) with the 
assumptions g^ (e )/y (e ) -* 0, e 3 /y(e) -• 0 and f 2 (e)/y(e)~0(l)ase-*0, we have: 


e^o7fe{ Ce31o8e) [T hx + ( 2 ^o-T h )] + f > (e) [ 2 ' J o Re(logZ) + T =] 


-g*(e)[cpi] } = 0 (VI-14) 

taking y(e) = e 3 loge, we have 

g 4 (e ) = e 3 log e (VI-15) 

f 2 (e) = e 3 log e (VI-16) 

g 3 (e) = log e f a (e) = e 3 log 2 e (VI-17) 


rcpil = -f (S)h (S)X+[2 U 3 (s)h (s) + T a (s)] 

L J o L 

+ 2 (ji 2 R e(log Z ) (VI-18) 
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(VI- 18) gives the inner boundary condition for CP 4 . We proceed to cp g . 
Again, we have the matching condition for m = 3, n = 5 as 


' SL Si cL 

^{ eS [(f h+P l) X+ T hXRe(1 °8 Z > + ( 2 U 3 -^ReClogZ) 
+ I 3 + -T-(I + |) + c 3 ]-8 s (e)[<p i B ] o }.0 


By taking Y(e)=e 8 , this gives 


g B (e) = e 2 


3. 3 

[cp s 3 o = (t ( s) h (s) +P 1 (s))x + T 3 (s)+G 3 (s) +[2\S 3 q - T l (s) h (s)]Re(logZ) 

a h(s) a x (s) . 

+ -f(s)h (s)XRe(logZ) + g (J +|J 

Now, the inner solution stands 

cp i ~cp^ + (e log e) cpi + e cpli + (e 2 log 2 s) cp* + (e 8 log e) cpl + e 2 cpt 


Substituting this back into the basic differential equation (III- 8 ) and 
the basic boundary condition (III-ll) give the governing differential equa- 
tion as 

V JW °P3 = 0 ( 


i dcpi 

V XY 'P 4 " h(s ^ bX 


i ^CPs r 3 s 3 -| i 

^cpb =h < s >-§r- iw 2 " h(s) _K 

and the inner boundary condition as 


V XY F _ i 


TvT x ^ 3 

\ 5n )c 

7 XY F _ i 



l ) 0 


(VI-27) 


31 


V XY F n i _r^\ _ hf/bs uc Po 

'vTlT * \bnJ “ |v vv f| ~5T 

AY a ay 1 


(VI-28) 


(VI- 28) can be rewritten as 

a 

where n is the unit normal vector to a in X-Y plane, dn represents 
a inf initesimile displacement in the direction of n and ^ = |-— |. 

c) Solutions for cpg and cp 4 

Since cp^ is independent of X, (VI-24) reduces to a two-dimensional Laplace 
Equation. Together with (VI-27) and (VI-18), cp 4 is governed by the following 


set of equations: 


V XY ^ = 0 


(VI-30) 


S cl 

— j(s)h(s) x + j2|j 3 - — -(s)h(s) + T S (s)J + 2^3 Re(logZ) (VI-30b) 

/Scp 4 x 

(sr) = 0 (Vi-300 


Compare this set of equations with the governing equations for cp* (V-28a,b,and c), 
we see that the solution for cp 4 is just the Re-entrant Flow Model with a local 
velocity at infinity as a^/2 instead of sina(s) for co \ . Actually, cp* is a 
correction to cp 2 » i.e., the local velocity at infinity, sin a(s), is only 
correct to 0(g) and we need a correction to this if we go to a higher order 
0(g 2 loge). Of course, if the jet trajectory were straight line [h(s) = 0], 
there will be no correction whatsoever. Thus, cp 4 represents part of three- 
dimensional effect of curvature on the inner solution. 

Define a complex w 4 (Z,s)as 

a Q (s)h(s) 

w 4 (Z.s) = 2 sina(s ) w 3 (Z,s) 
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(VI-35) 


a (s)h(s) co a 

o ~ • -zr T I sinaZ+a log Z +E r# 
2 sina(s) |_ o ° , Z n 


where w 2 (Z,s) is the Re-entry Flow Model defined in section (V-e). Then 
cp 4 and the source strength |_l® for c p° will be 


cp 4 (X, Y , s ) = Rej^w 4 (Z , s) - 


a (s)h(s) a 1 

Tl iT drG T b 0 (s) + ^o (s) - T <B)h(-)+ T ,(sa] (VI-36) 


a (s) h (s) 

^o (8) ■* 2 sin a (s) a 0 (s) 


(VI-37) 


cd 3 is governed by (VI-23, 26, 6). The solution of this set of equations 


is 


a (s)h(s) 

CP 3 (X»Y, s) - 2jj®(s) = 2 ° lno(s) »o< 3 > 

which represents another flow field tangential to the jet trajectory. 

d ) co^_ 

The differential governing equation (VI- 25) is Possion's equation with 
right-hand side to be a function of cp 2 and cp^ . The solution of Possion's 
equation can be written as a sum of homogeneous solution and a particular 
solution 

i _ i . i 
'Ps-cp 5h + 'P B p 

Let w 2 be the complex conjugate of w 2 , then cp °2 =i[w s (Z,s) +b Q (s) 

+ Tj^ (s) +w(Z) +b Q (s) +Tj (s)"j and (VI-25) becomes 

lift ■ ^ [*& + • h(s) slna(s) 

In this form, a particular solution co _ „ is easily obtained 

5 P 

cp g P (X,Y, s) = ^[zw a (Z)+Zw 8 (Z) - 2ZZsina(s)] 


(VI-38) 


(VI-39) 


(VI- 40) 


(VI- 41) 
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With the aid of power series expansion of w a (V-36), we have the behavior 
of cp p B at outer boundary of inner region as 

[tPB P ] - | [ Z b 0 + Z V »„(Z logZ +Z logZ) +«!(|+ !)" 


(VI-42) 


The homogeneous solution should satisfy following set of equations: 

s h ~ G 


= h ] „ - L<p 6 ] o - p] £ 

/Scpgh N = dn 5jP_o 
\ 3n / ds 3s \ 3n 


(VI-43a) 

(VI-43b) 

(VI-43c) 


Inserting (VI-21) and (VI-42) into (VI-43b) gives 

[cothl =(^ h + P i) X + T 3 +G 3+( 2 H® --^hjReClogZ) 
L - J o 

ah _\ 

+ -f- (Z logZ+Z logZ) -g (Zb o + ZbJ 


(VI-44) 


Attention should be called to the fact that up to now the source strength |j 3 
for cp° is still undetermined and is at our disposal. For convenience, cp g ^ 
and |j 3 are broken up into two parts: 


^ s h cOshi^Vsha 


2^_a 1 h/2=n 3 i +n 


3 

O 2 


(VI-44a) 

(VI-44b) 


where cp 5 ^ l i- s governed by 


^XY<P 5 hi = 0 


[cPs 1 hi] = [ c Psh] "Ho a R e(l°g z ) 


(VI-45a) 

(VI-45b) 


Here, the inner boundary condition for cp b h i ^ as not been specified. Of course, 
the solution to this set of equations will not be unique. By inspection, a 
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solution can be easily found as 


a (s)h(s)/2 + P 1 (s) 

cp Bh i(X,Y,s) = Re[ - ln - Q - (s) [w 2 (Z,s) -b o (s)]J + T 3 (s)+G 3 (s) 


a Q (s) h (s) 
8 


(Z logZ+Z logZ) -^^[Zb (s)+Zb (s)] (VI-46) 


e) 


i 

‘ •P 6 h a 

The governing equations for cp B ha are 


v ^Y<Pshs = 0 

L^Psha] = Re [^o a logZ ] 

dn^^o f 3 / i , i \~ 

= ds Ts" ' I^^bp +cp 6hl ; 
a o 

Let Yg ha t ^ ie stream function corresponding to cpg h 2 > or 


(VI-47a) 
(VI-4 7b) 


(VI-47c) 


w sha cPshs + *-^6ha (VI-48) 

Then we can obtain the solution for cp^ha hy solving a boundary value problem 
on Y 5 ha with the more convenient inner boundary condition 

a ^ i a 

Yghs (a,o) = -J d a +J[^(cpBP+cpghi)] (VI- 49 ) 

o o o 

where da represents a small displacement along the boundary a in Z-plane (Fig. 5). 
Due to complexity of the geometric shape of the boundary surface a» it is very 
difficult to solve this boundary value problem in Z-plane. Fortunately, there 
are conformal mapping techniques at our disposal and the surface a can be trans- 
formed into a geometrically simpler boundary. Through the shape in w s -plane 
and £-plane (Fig. 5) are much simpler than the shape of a in Z-plane, they are 
still somewhat complicated for our purpose. A better choice would be by the 
transformat ion 
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(VI-50) 


-Kc + z) 


to map the unit semi-circle in £-plane to the whole upper plane on r|-plane (Fig . 5), 
In ri-plane, the solution for cpghs * s eas ity obtained 


i i f ^ shs p -i 

Y S h 2 Cr 1r ,r li ) =- j ( vT)s+T1 jdT -H 0 3 a Inijlog (n- iC^) (n + iC^)] (VI-51a) 

- 00 

where T = dummy variable of integration 

r\ =ri r + iri i 

iC = the position of Point C in rrplane (Fig. 5). [Point C in Z-plane is 
11 

infinite.] 


The harmonic conjugate function of Y B hs : *- S: 

cpshs <‘n r ,n i ) - n 


1 f Y sha > °)^i1r ” T ^ 


(*n r -t) s +n? 


d T 


■ ^0 2 R e[ log ^ " iC n > (“H + ic „ ^ 


and 


00 1 


w Bhs <Tl) c P5ha +i '^5h2 n 


_ 1 i» ^ sha (tlj-’o) 


(r| -t) 


d T 


(VI-5 lb) 


U 3 g log (t, - iC n )( T1 + iC w ) 


il 


(VI-51c) 


By conformal mapping technique, w 5 h 8 (ri) can be transformed back into Z-plane 
to obtain w g h 2 (Z). However, we generally require not the potential but the 
velocity, i.e.: 

d Y B h a^ 

. .1 r-~^r (T,o) JT - u s r 1 - + L. - 

<3 Tj tt J (t]-t) 02 L (ti ~ iC ) (rj + iC ) 

-oo 1 il V 


^Bha ^ Bha d r| , dZ 
dZ dn dQ d £ 


(VI-5 2a) 


(VI-5 2b) 
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The next step is to evaluate d /dr| . Before doing this, we should 

clarify two more points. In section (V-d), we stated that mathematically the 
Re-entrant Flow Model consists of two sheets of a Riemann Surface with dotted 
lines EF and EF' (Fig. 5) as the boundary streamlines on the second sheet. 
Physically, we are only concerned with the flow field on the first sheet. 

All fluid flow across FF' is assummed to be carried away by the jet. But 
from a mathematical point of view, we have assummed that the inner boundary 
of our problem consists not only of the boundary lines ABHF and AB'H'F' on 
the first sheet but also of FE and F'E on the second sheet. More precisely, 
we have assummed - 1 

) = d7 TT "[^(vbP + cpBhi)] (VI-52c) 

a a 

where a denotes the boundary ABHFE and AB'H'F'E (Fig. 5). In essence, we have 
introduced artificial boundary condition on EF' and EF. One natural question 
will be: what is the effect of this boundary condition in the second sheet on 
the flow field of the first sheet. In order to answer this question, another 
boundary condition is introduced for comparison with equ. (VI-52c), i.e. 

— j =0 on EF and EF' (VI-53) 

a 

For Q=-1.5, the solutions to these two boundary conditions are calculated 
and found that there is essentially no difference between these two solutions, 
and therefore, the more convenient boundary condition (VI-53) is used in sub- 
sequent works. 

Another point to be clarified is how to determine . The potential 
(cp 5 P+cp 5 hi) will induce a flow field which is generally not tangential to 
the boundary a. cpghs i s introduced such that (cp B p + cp B ^ j+ep sh® ^ will be 


^<Psh2 
. S n 
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tangential to the boundary q . Therefore, cp 5 h a along will give a net influx 
of fluids through the boundary a- This flux must be supplied by "sources" at 
infinity in order that the equation of continuity is satisfied. But we have 
two infinity: one on the first sheet (C in Fig. 5) and another on the second 

sheet (E in Fig. 5) sources at either one of these can supply part of the 
flux across the boundary a. In order to determine \ 1 q 3 * following criterion 
is adopted: all flux across the portion of a on the first sheet is supplied 

by the infinity on the first sheet and all the flux across the portion of a 
on the second sheet is supplied by the infinity E on the second sheet, i.e. 


2nn 0 3 8 = -[V Bh2 (F,o) ■ Y Bhs (F 1 »o) 


(VI-53a) 


Now, we proceed to evaluate the boundary condition in rj-plane. Remembering 
S cp°/S S = cos a (s) , we can obtain the first integral of (VI-49) as 


ran 


SA. 


JS”s -aTda = -cosa(s) 77 


(VI-54) 


where is the shaded area in Fig. 6a. The second integral can be evaluated 

in the w 2 -plane. By using the analytic property of the mapping function w 2 (Z), 

it can be shown that 

a a 

• / > > \/ 4 -S \ 

(VI-55) 


J^^sP+CPshl)] d ^ = J i (§7: *5T 8 )( C PBP +C PBh 1 )<iW | 


T| 

U V 

Inserting cp B p (VI-41) and cp^i (VI-46), performing the differentiations and 
noting that w 2 = w 2 on the boundary u> we obtain 

* »h = fe.o) - - ^ Im{[[a o (l<,g Z + 1) & - Z ^ + 




-2Zsina^ -b o g ^dn}-cosa(s) -yf 


(VI-56) 
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II 


dC . . n 

" herc 7n - - 1 + ^trr ■ 

d.W n r\7 

— — - and — are defined in (V-31) and (V-34), (VI-45) can be integrated by 

d C d C 

by numerical method in either £ -plane or r] -plane. However, it is much easier 
to obtain dY^hs/drir numerically since no integration is involved in this case. 
To obtain dy ghs /d’r| r ('n r >o), we first rewrite the boundary condition (VI-49) as 


d^sha an^cp o r 3 / i i \*i 

~ fo,°) + 


(VI-57) 


Again, by using analytic properties of the mapping function w 2 ,we can prove 
i 


„ Qr l °*5h 2/ s „ .v( a d V i . i \ 
U a d^7 _ ^r ( ^ ,o)=U a 1 L ( ^ 2 ' ^ a A c PBP + cP6hx; 


cos a 


(s) 


(VI-58) 


after substituting cpgP and cpghi > this becomes (here is the velocity of 

Re-entrant Flow Model on ct) 

S Ysha 


4 In {f,[V losZ + 1) -V"* +Z ^E 




-2Z sina(s)]-^^^|S C o S a(s) (VI-59) 

This equation is evaluated by computer and is presented in Fig. 9 for Q = -1.5, 
-2.0, -3.0. 

From (VI- 43a), we have 

4 0 3 2 = -[ Y °6ha <F,o) - Y* ha (F',o)]/2n (VI-60) 

where Yghs I s defined by (VI-56) and points F,F' are shown in Fig. 4. If we 
examine (VI-56), we will see that only two terms are involved: one represents 

the effect of curvature h(s); another, the growth of the jet. Let 
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(VI-61) 


„ . cos a (s) 5a 

U 3 = a + — ~ v— 

p 0 2 OB 2 TT OS 


where A is the area enclosed by FAF'F (Fig. 4), and 


o 5 




(VI-61a) 


f) The source strength |j a (s) for cp 3 


Combining (VI-41), (VI-46), and (VI-51), we have the solution for co*. as 


cp 


a h 


5 = ^[ Zw 2 +Zw 3 _ 2ZZ sina-Zb o - Zb Q l + (Z logZ + Z logZ) 

r a 0 h / 2 + p i 

+ »«[ sinoT < , = -V + "sh = r I S +S > 


Appropriate outer expansion of cp 5 is 
,a h 


t . , d II . f . 

Cp 6 | = (-f- +P 1 )x + T 3 +G 3 + (m o 3 +M 0 3 J Re(log Z) 
~o 

+ ^XRe(log Z)+ i|3(| + |) 


where 


„a . (a h /2 + p , ) ■— r 


o 1 


sin a 


(VI-62) 


(VI-63) 


« 5a cos a 

^0 3 3 o 5 5s 2 tt 


Now, comparing (VI-63) with the outer boundary condition (VI-21) gives 


P, 


;<s) -[■ 


a x (s)h(s) a®(s)h(s) 

4 + 4 sina(s) 


a o5< s) n 


Pi (s)a Q (s) cos a (s) 5A 

ft • _ . / \ * » "V 


2 s in a ( s ) 


4tt 


5s 


(VI-64) 


This source reflects essentially three kinds of effect: the first term repre- 

sents correction due to curvature effect; the second represents the effect of 
the shape of jet trajectory, i.e., the addition to the local flow induced by 
the singularity distributions along the rest of the jet [function P l is defined by 
(V-5b)]; the third, the growth of jet cross-section. 
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-VII. COMPLETE POTENTIAL, VELOCITY, PRESSURE COEFFICIENT, 
ENTRAINMENT AND DRAG 

a) gp* - , q*~ and C^ 

Hie results of previous sections are assembled to give the complete 
inner solution up to 0(e 2 )as 


cp 


^•(X.Y.sje) = fcosa(T) dT + (e logs) a (s) 


+ £ Re[w 2 (Z, s) - b Q (s) + T x (s)] 


a (s)h(s) a (s) 
+ (e 2 log 3 e) 2 — 


sin a (s) 


f a (s) h (s) 

+ (e 3 log e) Re{ 2sIna(s , [w a (Z,s> - b o (s>] 


{ 


a 3 (s) h(s) 


2sina(s) os 


+ a (s) +— —r 


P x (s) a (s) 


sin a (s) 


5a cos a(s) 
3s 2rr 



+ e S {^8^ [z w 3 (Z , s ) + Z w 2 (Z , s ) - 2ZZ sina(s) 
-Zb o (s) -Zb o (s)j 
a (s)h(s) 

+ g (Z log Z +Z log Z) 

ra (s)h(s)/2+P 1 (s) 

+ Re L ? 3 in 0 (s) — [-.».■) -b„(.n +»,„,] 

+ T 3 (s)+G 3 (s)j- + — 
with. X,Y finite as e->o 


(VII-1) 
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where w 2 and w g h 8 are defined in sections (V-e) and (Vl-e) respectively. 
This complete potential for the inner solutions has not been evaluated 
because we generally require not the potentials but the velocities. To 
obtain the dimensionless velocity q*, we substitute the potentials in 
the following equation 



rfx B \ a , e s 
e e 1 - ehY 


3 1 i 


(VII-2) 


After simplf ication, we have 

rr a (s)h(s)-, - n 

q^x.v.s-.c) ~ e x {[l+(e lo 8 e) 2 s i na ( s ) j u 2 + e U B 1} 

r a (s)h(s) 

+ e Y {j_l+ (e log e) 2sin “ (s) Jv 2 +e V B 

r da pScp 2 

+ e s |cos a(s) + (g log e) (s) + e | +h(s)X cos a(s) 



+ ... 


with X, Y finite as co 


where 


dw. 


U 3 - iV 2 = (Z , s) , complex velocity of Re-entrant Model 

Gw 

.. _ _ _ dw 2 

U 5 - i V 5 = — ^ S ^ | a 0 (s)[ log Z + ll-b o (s)+w 2 (Z,s)+Z -^-(Z, s) 

2 a (s) dw g 
- 2Zsina(s) + s . na(s) 

P x (s)a Q (s) 

+ sin a (s) 


dw 


5 h 2 


dZ 


(Z,s) 


(VII-3) 


(VII-3a) 


42 



From (VII-3), it is seen that q* consists of four elements: 1) the basic 
two dimensional solution dw 3 /dZ; 2) corrections due to the three dimen- 
sional effect of local curvature, which is reflected by all terms con- 
taining h(s) and some part of dWg^/dZ; 3) corrections due to the growth 
of jet cross-section, which are contained in the term dw 5 ^ a /dZ; 4) Higher 
order corrections due to the singularities distributed along the trajec- 
tory, which are represented by the function P x (s). 

The pressure coefficient will be 

Cp(X,Y, s:e) = 1 -q 1 • q 1 

= [l-(Uf+V§)- cos 2 a(s)] 

a (s)h(s) da 

' 2 <e lo © e >[ 2 sinG(V) <"! ) + cos a (s) -jj- (.)_ 

r fa 8 \~i 

■ 2eJ^(U a U s + V s V g )+ cos Q: ( s ) (^-^- + h(s) X cosa(s) j 

+ . . . 

with X,Y, finite as e-*o (VII-4) 


, . o , _o 
b) cp and q 

The complete outer solution up to 0(e 2 ) is: 

s 

cp°(x,y,s:e)~f" x sina(s) +J* cos (t) dTj 


s 


a (t )/2 


+ e 


K, 1 


d T 


+ (e ® log e) J 


a 2 (t )h(T )/4 sin a(T ) 


d t 
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with x,y fixed as e 


(VII-6 ) 


where 


M-(S) = a Q (S)/2 + ( e logs) a » (S)h(S)/4 sina(S) + ^ (S) 


is defined by (VI-64) 

It is seen that q° is the superposition of three elements: 1) the basic free- 
stream cross-flow Ug, i ; 2) a line distribution of doublets representing the 
local two dimensional blockage effect; 3) a line distribution of sources: 
the first term of the source strength reflects the effect of entrainment. 

The second and third term represent correction due to effect of local curva- 
ture, the growth of jet cross-section and the shape of jet trajectory as 
explained in section (Vl-f). 
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c) Total Entrainment E(s) 


The mass entrainment E(s) of the cross-flow at each station s is 

assumed to be the two dimensional flow-rate of the re-entrant flow, which 

can be expressed as / 

r 

E(s) = - f q 1 ’ n do (VII-7) 

v 

F 

The curve FT' are shown in Fig. 4. Examining q* defined by (VII-3) and 
recognizing 

f' 

- I" (U 3 ®x + V s ®y ) • n do = Eg (s) , 

F 

we have 

a (s)h(s) 

E(s; e ) ~ E 2 (s)+ (e log e) 2 "sina(s 7 E s( s ) + eE s( s ) + 


where 


with X,Y finite as c-*o 


(VII-8) 


E 2 (s) is given by (V-50) 


r t p i _ a (s)h(s)/2+P x (s) 

E s (s) = -| (V x ^cp 5 ) *n do = -| (’ XY tp 5p ) ' n do + ^ n " a (s) E s (s) 

F F 


Since E 2 (s)/R.(a) sina(s) is a function of cavity number Q only, E(s)/Rj (s)sin(s) 
is a function of cavity number only up to 0(1). We have assumed Q to be con- 
stant along the trajectory; therefore, E(s) will be proportional to 
Rj (s)sin a(s) to 0(1). This means that E(s) will have a maximum near the 
jet orifice and diminish gradually as the jet trajectory become parallel to 
the free cross-flow. This corresponds qualitatively at least to the physical 
situation. 
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d) Total Drag D(s) 


The drag D(s) is determined by using momentum principle (Ref. [16]) 
applying to the system shown in Fig. 10. AB and A / B / are two infinitely 
long plates symmetrically located with respect to X-axis. Since sections 
AA' and BB' are assumed to be located infinitely far away upstream and down- 
stream, all disturbance created near the jet are subsided and the velocities 
q. and q„ across these two sections are uniform. 

By applying continuity equation to this system bounded by AA'bb' and a, 
we have 

P E(s) = (q A -q B )pi (VII-9) 

and the momentum equation gives 

F' 

(q B " q A } pjJ "I (pq X )( %Y ‘ “> da = (P A " P B )je " D (VII- 10) 

F 

where q x is defined by the equation 

5 1 ■ Vx + Vy + V. 

q XY = q x e X + q Y e Y 

Combining (VII-9) and (VII-10), and using Bernouli's equation 
(P A" P B )i= ^ p(q B S “ q A 3) 

we have 

F' 

D = | p E (q B + q A ) + J (P q x ) q x dy (VII-11) 

F 

Now, let £-co, then q B -sina(s), q A -sina(s), we have 
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(VII-12) 


- pE(s) sina(s) + f p q| dy 


Substituting E(s) (VII-8) and % , Vs+ < £ Iog £) 
into this expression gives; 


a 0 (s)h(s) 

2 s in a ^2 + £ 


0(s;e)~p sina(s) e 8 ( s ) + ( e iog e ) . 

2 sin a;(s) Sln a ( s ) P E s ( s ) + epsin a(s)E s ( s ) 
+ .f[ U § + 2 ( g logs) ~ .- )h(S) - U2 +2cU a *] 

p L 2 sm a(s) 2 ^ cu s u s j p dy (VXI-13 

on the curve FF', is approximately equa! to (. „ > „ hI . 

4 ^ u w h which is a constant, 

and recognizing that 

F' 

J U e dy = "E g ( s ) 


D(s;e) ~p E 2 (s)si n a(s) +u 1 + r, Ioo.vp * 0 < s)h < s > r 

L vj lel0ge>()E o( 8 >27I?H(I)L Slna(s) + 2 > I „ 

+ epE 5 (s) sin a(s) + 2 uj + . . . 


with X, Y finite as 


(VII-14) 
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VIII. JET TRAJECTORY 


a) Cavity Number and Trajectory 

Up to now, a very important parameter remains unspecified, i.e., the 
cavity number Q. We have seen that the geometrical shape of inner boundary, 
entrainment and drag are all functions of Q. If we can derive a relation- 
ship between these three elements, then Q will be determined, and, in turn, 
geometrical shape, entrainment and drag will be determined accordingly. 

The trajectory will serve to give this relationship. 

The curvature of the trajectory is directly influenced by the combined 
effects of pressure and shear stresses acting at the jet boundary. These 
stresses not only change the momentum of the original jet fluid but alter the 
momentum of a portion of the free stream fluid and through the agency of 
turbulent mixing give rise to entrainment as well. In the inviscid model 
employed here the role of viscous stresses has been reduced to that of 
accounting for the origions of the vortex sheet a w and the resulting wake 
region while turbulent mixing has been replaced by the re-entrant flow. One 
of the principle advantages of the model employed however is that it provides 
a dynamic coupling between the entainment and the pressure stresses typified 
by the cavity number Q. As a result, the value of Q is a determining factor 
in the trajectory shape, for which Margason (Ref. jlO}) has provided experimental 
data. Thus, trajectories are computed for several values of Q and compared 
with Margason' s trajectories. 

The value of Q yielding the closest fit is then chosen for the computa- 
tion of pressure coefficients. 

b) Derivation of the Trajectory 

In section (V-d), we assumed the front surface of the jet (cij, Fig. 4) to 
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be a circular arc with radius Rj(s), and rear surface BFF'B' unspecified. 
Here, we still let the shape of the rear surface remain unspecified, but 
the total cross-sectional area is assumed to equal an equivalent circle, 
e. i. , ttR? (S) . 

For an inf initesemally thin cross-sectional element (Fig. 11), the 
momentum equation and continuity equation can be written as 


dM j 
dS 


D(S ) e x + pE(S) cosa(s) e g 


(VUI-la) 


dS = pE(s) 


(VUI-lb) 


where m^(S): The rate of mass flux at station S 


m . = ttR f P v . 
J J J 


(VIII-2a) 


Mj(S): The rate of momentum flux at station S 


M. - m.v.e s .nE?pv.e g 


(VIII- 2b) 


Vj(S): Velocity of the jet at station S 

S: Dimensionless coordinate along trajectory S = s 1 /R^ o 

Here, we have taken the average velocity of the jet as Vj(S)(see Fig. 11). 

Substituting M. and m. into (VIII-1) and simplifying it, we arrive at: 

J J 


da _ -D(S) 

dS ttR? (S) Vj(S)p 


(VIII-3a) 


dV j _ E(S) f cos a (S) 

"3T ttR^CS) L Vj (S) J 


(VIII-3b) 


dR j = E (S) r 2 cos a (S)- [ 

dS 2 ttVj(S)Rj(s) L v j (S) J 


(VIII-3c) 
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Now, substituting E (VII-8) and D (VII-14) and after some algebraic work, 
we finally have 


. E sin 2 a(S) , 

da _ _o r ah E B 

dS _ frRj(S)v 2 (S) L(l + U o ) + (e loge) 2 ^ (1 + 2U q ) + e (1 + 2U q )J (VIII-4a) 


dv. E sina(S) . . s _ 

1 o ,'cos a(S) , > f i , / , 

^ ' -^OsT - - i/[i + <€io ge ) 


a h 
o 


2s in a 


+ e 


-1 

E = J 


dR . E sin(S) 
J _ _o 


dS 


2rr v 


n(S) cosa(S)\r . a ° h E s-| 

71 s) \ 2 --TasT) L 1 (e loge) 2 + e TL J 


with the initial conditions 


where 


a( o) = tt/2 
Vj(o) = IT 

Rj (o) = 1 


E = 
o 


R . s in a 
J 


= TTb 


d 3 1 

74 a-Q) 2 


a function of Q only, see 

U 


U 

o 


w 

sin a 


(V-50) 

(1 -Q)* 


(VIII-4b) 


(VIII-4c) 


(VIII-4d) 

(VIII-4e) 

(VIII-4f) 


also a function of Q only. 

This is a set of coupled nonlinear first order differential equations. We 
cannot integrate them immediately, because we do not know e beforehand. But 
we can solve it by iteration methods. First, take only terms up to 0(1) 
i.e., the first term in the brackets, integrate them by suitable numerical 
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methods to get a(S), R^(S) and Vj(S). Then we obtain the actual trajectory by 
integrating numerically the following set of differential equations 


dX 

= cos a(S) 
dZ 

~dS = Sina(S) 

with the initial conditions 

X (o) = 0 
o 

Z (o) = 0 
o 


(VIII-5a) 


(VIII-5b) 


(VIII-5c) 

(VIII-5d) 


where X , Z are dimensionless Cartesian coordinates, see Fig. 11. 
o o 

From the solution of this set, e> h(S), PjCSJare obtained and substituted 
back into (VIII-4). Now, we can take the full equations and integrate them 
again. This process is continued until the difference between two successive 
solutions are small. For the value of U^/U^ considered here, the vary first 
integration was accurate enough for our purpose. In Fig. 12, trajectory for 
Uos/U^ = .8, .5 and .25 are plotted with Margason's curve. The jet velocity 
Vj(S) and jet radius Rj(S) along the trajectory are shown in Fig. 13 and 
Fig. 14, respectively. 


51 



IX. CONCLUSIONS 


The inner and outer solution for pressure coefficient for Uoo/U.. = .8, 

.5 and .25 are presented in Fig. 18, Fig. 19 and Fig. 20, respectively. 
Experiment data by Bradbury and Wood (Ref. [7]) and by Vogler (Ref. [8]) 
are also presented with the outer solutions for comparison. 

Since the reentrant flow model introduces a stagnation point behind 
the jet, it is incapable of reproducing the experimental pressure contours 
in the wake region. However, upstream of and along the sides of the jet, 
the computed pressure contour gives good agreement with the experimental 
data for Ufco/Uj = .25. For U m /Uj = *5 and .8, the computed pressure contour 
seems to have the right order of magnitude but is not swept back enough. 

It is possible that the inclusion of higher order terms would remedy this 
as indicated by the comparison of first order and second order solutions 
in Ref. [4] . 

As indicated in the introduction, Werner and Chang (Ref, [4]) raised the 
possibility that representing the effects of entrainment by a net sink dis- 
tribution along the jet trajectory might not be correct. In order to 
answer this question, we have calculated the contributions to the source 
strength from the various elements of the flow. In our solution, the 
outer flow region is represented by a line distribution of doublets and 
sources. The doublets simulate the effect of "local blockage", while the 
sources, as we have repeatedly stated, are derived from four different 
origins: 1) the entrainment, 2) the growth of the jet cross-sectional 

area, 3) the curvature effect, and 4) higher order correction due to 
singularities distributed along the jet trajectory, i.e., the effect of 


Pj (s) (V-5b). The magnitude of these contributions as a function of Sj/R^ 
are shown in Fig. 15, Fig. 16 and Fig. 17 for iWUj = .8, .5 and .25, re- 
spectively. From these results, the following conclusions are drawn: 
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% 


1) The source strength due to PjCs) has been calculated and proved 
to be very small for all cases. Actually, it is too small to be shown 
on these graphs. 

2) Around Uj/TJco = .8, the combined effect of entrainment and jet 
growth is small. In other words, if only the source due to curvature is 
used, a good approximation is obtained. 

3) At iWUj = .25, the effect of curvature is much smaller than the 
combined effect of entrainment and growth of the jet. 

Judging from the trend, we can neglect curvature effect for He/lL £ .125. 


4) Introducing 




(IX-1) 


we see that the source strength due to entrainment and the source strength 
due to growth of the jet tend to cancel each other for 

•!%><•!/ Vo 

5) It seems that, for s 1 /R^ o < (s x /R ^ q ) o , following two equations 
will give a reasonable approximation 

- f ° r — s 4 - 5 Cut) - 3 <“- 2 > 


H e = Cl 


J° 


s i s i / U i \ 

!t *r ' £or — s 4 - 5 lui )' 3 

6 JO JO 


(IX-3) 


where C x , C 3 are constants 

n : source strength due to entrainment 

^e 

: source strength due to growth of the jet 

8 
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Fig. 1 - Jet In A Cross-flow And Mathematical Model 




JET TRAJECTORY 





Fig. 3 - Local Coordinate System 







Fig. 4 - Re-entrant Flow Model 
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Fig. 6a - Contour a 
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Fig. 6b - Contour <j> Q = -2.0 
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Fig. 6c - Contour <j, Q = -3.0 



Fig. 7 - Entrainment E 2 and Drag Coefficient 
vs. Q, Re-entrant Flow Model 




Fig. 8 - Position Of Rear Stagnation Point vs. Q, Re-entrant Flow Model 







Fig. 10 - Control Volume For Determination Of Drag D(s) 
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Fig. 11 - Derivation Of Jet Trajectory 




Fig. 12 - Jet Trajectory 






Fig. 13 - Jet Velocity V 





Fig. 14 - Jet Radius 
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Fig. 15 - Comparison Of Source Strength Of Different Origin, Un/U 
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Fig. 17 - Comparison Of Source Strength Of Different Origin, Uoo/Uj = .25 






Fig. 18a - Pressure Contour, Outer Solution, U ro /U . 



Fig. 18b - Pressure Contour, Inner Solution, U /U< 
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Fig. 19a - Pressure Contour, Outer Solution, Cto/U. 
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Fig. 20a - Pressure Contour, Outer Solution, U^/U 
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